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1 . Introduction 

Statistical mechanics provides a molecular founda- 
tion for thermodynamics. This results in the expres- 
sion of thermodynamic properties as functionals of 
the intermolecular potential functions of the constit- 
uent molecules. In principle, given the potential 
function appropriate to a given system one can cal- 
culate all of the thermodynamic properties of that 
system merely by turning the computational crank. 
In practice, matters are not so simple both because the 
relationships to be evaluated are enormously compli- 
cated and because the potential functions are not 
known with sufficient accuracy. 

The most accurate of the statistical mechanical 
expressions contain /V-body potential functions 
(7V—10 23 ) which are impossibly difficult to calculate. 
Simplifying assumptions can be made which often, 
as in the virial expansion, result in a power series in 
some parameter (e.g., the density) whose coefficients 
depend on lower order /V-body potentials (/V=2, 
3, . . .). Since even the three-body potential is ex- 
ceedingly difficult to calculate, the assumption of 
pairwise potential additivity must generally be included. 
With this assumption, the thermodynamic properties, 
in the statistical mechanical expressions, can be made 
to depend on only the pair potential function. Most 
theories for the further simplification of the compli- 
cated expressions proceed from this point on the as- 
sumption that the pair potential function is known. 
These theories result in simpler, but generally still 
complicated, relationships between the thermody- 



namic properties and the potential function. The 
evaluation of such theories can be a formidable task. 
All too often their final evaluation is obscured [l] 1 
considerably by the fact that pair potential functions 
are, in practice, only imperfectly known. 

The intermolecular potential function together 
with a sufficiently accurate theory can be used to 
extrapolate far beyond the bounds of available experi- 
mental data [2], something which is not possible using 
completely empirical methods. Such extrapolations 
are very strongly dependent on the potential function 
and can be considerably in error when the wrong po- 
tential function is used. 

Clearly the determination of accurate intermolecu- 
lar potential functions is of some importance. In 
this work we have sought to understand some of the 
methods generally used to determine such functions, 
particularly with respect to the question of the unique- 
ness of the potentials obtained. A lack of uniqueness 
exists when a set of experimental data for a given 
property can be correlated equally well using the ap- 
propriate theory and any of two or more potential 
functions. Where a lack of uniqueness exists, it 
becomes necessary to attempt to determine if there 
is a lack of sensitivity inherent in the theoretical 
quantity itself or if it is lack of experimental precision 
which makes it possible to fit the data equally well 
with two or more theoretical curves. These are 
equivalent to the following two questions. First, 
how well can the property of interest, in principle, 
distinguish among potential functions and, second, 
how well can it distinguish among such functions at 
the present time, given present day experimental 
precision. The answer to the first question is a perma- 
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nent one while the answer to the second one changes 
as experimental technique is refined, approaching 
the first answer in the limit of zero experimental 
error. We have restricted ourselves to the first 
question and discuss the second one only in passing 
mainly to place our results in a practical light. The 
answer to the first question is of considerable practical 
importance since it points out where refined experi- 
mental techniques will not produce more information 
about the potential function. 

In principle the pair potential functions can be cal- 
culated in an a priori fashion using quantum mechan- 
ics by calculating the potential energy of two molecules 
as a function of nuclear spearation. Since one needs 
to consider all the electrons in each molecule, this 
is also an intractable TV-body problem, TV now being 
the total number of electrons involved. This leaves 
one no alternative but to turn the problem around 
and determine the potential, in some manner, from 
experiment. In practice, the procedure is reduced 
to a semiempirical one. A functional form is assumed 
for the potential whose choice is based, in part, on 
theoretical arguments. In this form are included 
parameters whose values (and hence the detailed 
potential) are to be determined from experiment. 
The parameter determination is made by substi- 
tuting the potential into statistical mechanical expres- 
sions for some macroscopic property and comparing 
the result with experiment [3]. Best results are to 
be expected when the theory is one whose depend- 
ence on the pair potential is strongly based as is the 
case, for example, in the low density limit of certain 
theories. 

In this paper we shall discuss the use of the zero- 
density viscosity, diffusion, and adiabatic Joule- 
Thomson coefficients and the second virial coefficient 
as ways of obtaining the potential parameters. For 
each of these, the pair potential appears in the inte- 
grand of an expression for the macroscopic property. 
This suggests that, given enough experimental data 
of sufficient accuracy, one might be able to invert 
the theoretical expressions and determine the po- 
tential as a unique functional of the experimental 
data. For the second virial coefficient, however, Le 
Fevre and Keller and Zumino [4] have shown that the 
potential is not determined uniquely by the data, even 
in principle. In this work we demonstrate this lack 
of uniqueness quantitatively for particular potentials. 
We show there exists a lack of uniqueness for each 
of the macroscopic properties considered, although 
it is somewhat less pronounced for the transport 
properties than for the equilibrium properties. For 
each property, the lack of uniqueness is found to be 
more pronounced in one temperature range than in 
another. Strong positive statements can then be 
made about the temperatures at which experiments 
designed largely to determine potential functions 
should not be performed for particular substances. 
We are also able to show the simultaneous fit of cer- 
tain of these properties to be sensitive to differences 
in the potential function. 



2. Computational Method 

Because of its relative simplicity, we shall use the 
second virial coefficient to illustrate the details of 
the computation. The other properties are handled 
in essentially the same manner. We shall restrict 
the discussion to two parameter potentials. This is 
no real restriction since a three-parameter potential 
can be treated as a family of two parameter poten- 
tials one for each value of the third parameter. 

The second virial coefficient is related to the poten- 
tial function, <£>, by the relation [3] 



B(T) = b B*{T*)-- 



% exp 

J o . 






"* 2 dr 



(1) 
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where b = — , T* = kT/e, <p* = <ple, and r*=r/cr. 

Here, as usual, a is a characteristic length related to 
(p and 6 is the depth of the potential well. TV is 
Avogadro's number. For a given potential, a given 
7", and in the left-hand side, a given experimental value 
of B at that 7 1 , (1) contains only the two unknowns a and 
e. In this work, in place of experimental B(T) values, 
we supply to the left-hand member of (1) the second 
virial coefficient for a potential function other than the 
one appearing in the right-hand side. Thus, if the 
subscripts 1 and 2 refer to the two potentials, (1) 
becomes 



B(T) = (b ) 1 Bf(Tf) = (bo)2B*(T: 



or 
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exp 



T* 



- 1 \r* 2 dr* 



where B ?(T?) is a functional of <pf (r*). Now T} = kT/ei 
and T* = kT/e-i hence, for a given T, T* =^T* so that 



■l\r* 2 dr*. (2) 



This is an equation containing two unknowns, namely, 
the ratios 62/61 and (60)2/(60)1- These ratios are com- 
pletely determined, given another equation connecting 
some other functional of the potential for the two poten- 
tials. We have taken, for this second functional, the 
first derivative of B with respect to the logarithm of T 
[5]. We thus determine, at each temperature, that 
pair of ratios, e 2 /6i and (60)2/(60)1 which results in the 
equality, for the two potentials, of both the second 
virial coefficient and its first derivative with respect to 
the logarithm of the temperature. Ordinarily, the 
ratios are different at different temperatures, hence 
we shall actually determine their temperature depend- 
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ence. By dividing the equation which represents the 
requirement of equality of the first derivatives by that 
representing the equality of the second virial coeffi- 
cients themselves, one obtains a single equation for 
the equality of the logarithmic derivatives. The work- 
ing equation becomes 



S FOR B, T < T B 



d]nBf(Tt)_d]nB*(T$) 



dlnT* 



d\nT* 



(3) 



This is a single implicit equation in the ratio e 2 /ei, the 
dependence on (60)2/(60)1 being removed in the process 
of differentiating the logarithm. 

On the substitution of the set of values Bf{Tf) in 
the left-hand member, (2) becomes a family of curves 
for the ratio e 2 /ei, as a function of the ratio (60)2/(60)1, 
one curve for each value of Tf [7]. These can also 
be used for obtaining information about the potential 
function [8], although the results so obtained cannot 
be easily presented in a manner suitable for our pur- 
poses. We have therefore used a different approach. 
We have computed, for each of a number of potentials, 

fi* and B'* = r* (4^:1 and from these 



s= 



IT 

T* dB* = B'* = d In B* 
B*dT* B* din T* 



These are used to solve (3) in the following way. Each 
value of To has associated with it a value S 2 . Using 
inverse interpolation, the value of T* is found for which 
Si=S 2 . The ratios e 2 /ei and (60)2/(60)1 are then com- 
puted simply from 



rri ^ ano 
€1 T$ 



(60)2 
(60)1 



B$(T$) 

Bf(Tf)' 



This procedure can be illustrated graphically with the 
help of figures 1 and 2. These contain plots of S 
for the second virial coefficient [9] versus the loga- 
rithm of the reduced temperature for several rep- 
resentative potential functions. Temperatures below 
the Boyle temperature appear in figure 1, those above 
that temperature in figure 2. Note that the S values 
for each potential are plotted against the logarithm 
of the reduced temperature for that potential. Equa- 
tion (3) is solved graphically at each Tf by measuring 
the horizontal distance between the ordinate associ- 
ated with Tf on curve 1 and the same ordinate on 
curve 2. This distance is then just 

log 77 -log 77 = log (71/77) = log (€l/€ 2 ). 

It should be noted that where there are values of S 2 
greater than the maximum Si value, solution of (3) is 
impossible. This occurs (see fig. 1) for the (9,6) poten- 
tial compared to the (12,6). 

Note that both the value and slope of Bf are fit by B£ 
at Tf. It follows, therefore, that Bf will actually pro- 
duce a relatively good fit to B f in a small neighborhood 




FIGURE 1. S for the second virial coefficient of several potentials, 
T* less than the Boyle temperature. 
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Figure 2. S for the second virial coefficient of several potentials, 
T* greater than the Boyle temperature. 



around Tf for the parameter ratio obtained by solving 
(3) at Tf. Furthermore, when this ratio is found to be 
independent of Tf over a range of values of Tf, B* 
provides an excellent fit to Bf over that range. 

In this calculation, the Boyle temperature serves as 
a natural dividing point between two temperature 
ranges. Since B is zero at that temperature, the deriv- 
atives in (3) are not defined. Furthermore, since B 

changes sign while -pp does not, S also changes sign. 

Hence, values of S taken from temperatures above the 
Boyle temperature for one potential cannot possibly be 
made to fit those taken from below it for the other 
potential. Therefore, we have treated these two 
ranges separately coming as close to the Boyle tempera- 
ture as desired from either side. 

As mentioned in the introduction, the zero density 
adiabatic Joule-Thomson, viscosity, and diffusion 
coefficients are also considered here. These are also 
treated as outlined above. The equation (3) for each 
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of these is replaced in turn by the requirement that the 
quantities ^ for the adiabatic Joule-Thomson 



coefficient, 

d\n[&T* 
dlnT* 



-3/2 



dlnrjT*- 1 ' 2 
dlnT* 



shown that 

dlnT* 



=(/+2) 






-C/+2). 



for the viscosity, 



md 



for the diffusion coefficient be equal for 



the two potentials. Now [3] /x*C^ ) = T 



dB* 

dT* 



-B 



= B'* — B* so that the connection between this quan- 
tity and the potential function is essentially the same 
as between the second virial coefficient and the poten- 
tial. From this last, 



T*-^[vlCW]=:B"* 



d 2 B 



and 



d\n[fiC^] = B' 
d\nT* 



From this it follows that, for the viscosity, eq (3) is 

(|(2,3)* 

replaced by the requirement that 4 n(2 2) ^ — 4 be equal 
for the two potentials. For the diffusion coefficient, 

£}(1,2)* 

on the other hand, the quantity 3 n(1 1)s|c — 3 must be 

equal for the potentials. Since J) and 3) are each 
nonzero, there is no dividing temperature analogous 
to the Boyle temperature for these properties. The 
S values corresponding to these properties are plotted 
in figures 3 and 4. 



B 1 



■B* 



where B"* = T* 2 -^r-. Here, as with the second 
off* 2 

virial coefficient, there is a temperature which divides 

the calculation into two parts. Now, however, the 

zero density inversion temperature (defined by 

#'* =B*) plays the role previously played by the Boyle 

temperature. 

The connection between the viscosity and diffusion 
coefficients and the intermolecular potential func- 
tion is contained in the so-called collision integrals 
[3] 

n«>*)*(r*) 



(S + l)!T* (s+2) 



J n 



e -4**IT* g*(28+3)Q(l)*(g*)dg* 



where 

Q m (g*)= — ii 2 - f (-i) f / o aC(1 ~ cosi x)^*^* 



i 



2 1 + / 



with the intermolecular potential function being con- 
tained in the equation for the scattering angle 



X(g*,b*)- 



TT~2b* t 

J 1 



dr*/r* 



Vl-6* 2 /r* 2 -(/)(r*)/g-* 2 



r„\ being the distance between a pair of molecules at 
the time of closest approach. 

In terms of these, the zero density viscosity and 
diffusion coefficients can be written [3] 



^7^-1/2: 



C a 



Q(2,2)*(J*) .2 



3)T*-^ = 



C b 



jj(i,D*( t*) 



where C a and Cb are constants whose precise values 
are of no particular interest here. It can be easily 




LOG T* 

FIGURE 3. S for the viscosity coefficient of several potentials. 
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FIGURE 4. S for the diffusion coefficient of several potentials. 
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The macroscopic properties of most experimental The Square Well: 
systems have been correlated using the (12,6) poten- 
tial function [10]. As a result e/k values for that po- 
tential can be found in the literature for just about 
all possible systems of interest. Furthermore, it 
has been shown that, for very many substances, a 
very simple relationship exists between the e/A value 
for the (12,6) potential and the temperatures for the 
critical and normal boiling points for that substance where r* = r/cr. 
[3, 10, 11]. As a result, good first guesses for the 
efk values for the (12,6) potential can be made for 
any system for which either the critical temperature 
or normal boiling point is known. For these reasons 
we have chosen to compare each potential with the 
(12,6) function using the reduced temperature for the 
latter as the reference temperature Tf. The conver- 
sion to real experimental temperatures for any system 
merely requires multiplication by the (generally avail- 
able) e/A value of the (12,6) potential for that system. 



?*(r*) = 
<p*(r*) = 



00 f 



1, 1 



R* 



0. 



R* 



(4) 



3. Results 

The Second Virial Coefficient. Our results for the 
second virial coefficient are given in figures 5 to 8 
as plots of the ratios 62/61 versus the Lennard-Jones 
(12,6) reduced temperature. Figure 11 contains plots 
of the ratio (60)2/(60)1. The potential functions con- 
sidered are the following: 
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FIGURE 5. The parameter ratios 62/61 for the second virial coeffi- 
cients of the (m,6) and square well potentials with respect to the 
(12,6) potential. 

Note particularly the flatness of the curves for the (m,6) potentials. 



where r* = r/cr, a is that value of r for which <p* — 0. 
The Kihara: 
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- and 7* = — , a being the core 
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Figure 6. The parameter ratios 62/61 for both the second virial 



where r* = rlr m , r m being that value of r for which a ? d viscosity coefficients of the exp-6 potential with respect to 
# , the (12,6) potential. 

if 1 • Note that the former are flat while the latter are not flat. 
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FIGURE 7. The parameter ratios €2/61 for all four properties for 
the Kihara potential with respect to the (12,6) potential. 

Note the absence of any single ratio for which both the equilibrium and transport prop- 
erties are flat. 



Figure 5 contains results for the (mfi) and square 
well potentials, while results for the exp-6 and Kihara 
potentials are contained in figures 6 and 7. Results 
for the (12, ra) appear in figure 8. Each curve repre- 
sents the comparison between a potential and the (12,6) 
potential, the subscript 2 always referring to the for- 
mer potential. For each curve there is a smooth 
transition through the Boyle temperature indicating 
that the singularity in S at that temperature presents 
no problem. There exists, for each potential, a tem- 
perature range around the Boyle point such that in that 
range e 2 /€\ is essentially independent of 7 1 *. This 



means that, for temperatures in that range, the (12,6) 
and the potential with which it is compared are equiva- 
lent [12]. In fact, each curve is flat to such an extent 
in this range that a choice cannot be made between 
the (12,6) potential and the one with which it is com- 
pared based on experimental second virial coefficient 
data taken entirely within the range, even when these 
data are obtained with an impossibly high precision. 
What is particularly striking is that there is a single 
temperature range in which all the curves are flat. 
This range becomes exceedingly large if one does not 
include the square well potential in the comparisons. 
Obviously, there is a reduced temperature regime in 
which the second virial coefficient is particularly use- 
less as a probe of the potential function. What these 
results show specifically is that the second virial co- 
efficient cannot be used in this range to distinguish 
among any of the members of the (m,6), exp-6, Kihara, 
(12, n) and square well families of potentials. The 
list would presumably have been broadened had we 
considered other classes of functions. 

Of considerable interest are the results obtained 
when the (12, n) potentials are compared with the (12,6). 
The curves obtained for these are essentially the same 
as those for the (mfi) emphasizing the fact that the 
second virial coefficient cannot be used to determine 
the exponent of the attractive part of the potential. 
The requirement that the attractive exponent be 6 is, 
rather, a restriction placed on the potential based on 
a priori information, at least for the second virial 
coefficient. Thus, we see that the second virial 
coefficient is determined by the general shape of the 
potential and not necessarily by its details. This 
has previously been demonstrated formally by Le 
Fevre [4] and by Keller and Zumino [4J. They showed 
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Figure 8. The parameter ratios e-z/ei for both the second virial 
and Joule-Thomson coefficients for the (12,n) potential with re- 
spect to the (12,6) potential. 

Note that for temperatures above the Boyle point all curves are flat making it impossible 
to select a value of n at such temperatures. 
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that all potentials for which the sum of a certain pair 
of integrals, one over the repulsive part and one over 
the attractive part, were equal yielded the same second 
virial coefficient. A special case of their result is 
that all potentials with the same repulsive part and 
whose attractive parts have the same width as a func- 
tion of depth (i.e., but whose bowls are possibly dis- 
placed laterally) yield the same second virial coefficient. 

At temperatures outside the flat portion, the ratio 
62/61 is no longer constant. Below T : = 2.0 a par- 
ticularly rapid variation is indicated. The origin of 
this rapid variation can be seen in figure 1 to be due to 
large differences in the slopes of the corresponding 
S curves for equal values of the ordinate. According 
to figure 1, there are also potentials for which there 
is a temperature range in which there are S values 
larger than the maximum S value of the (12, 6) potential. 
The (9,6) function is an example. In such cases, 
solution of (3) is impossible. In other words, no ratio 
62/61 exists by means of which one can obtain simul- 
taneous equality of both B and T -7™ for the two 
potentials. 

The existence of a rapid variation of 62/61 with T ] 
at low temperatures would seem to indicate a very 
strong sensitivity, at such temperatures, to differences 
in the potential functions. Inability to solve (3) 
indicates an even stronger sensitivity to such dif- 
ferences. However, the sensitivity indicated applies 
strictly to exact data and the exact simultaneous fit 

of B and T — • As the requirement on the exactness 

of the fit is relaxed, the sharpness of the variation of 
62/^1 with F* is reduced. The introduction of these 
uncertainties in effect replaces each 5 curve of figure 
1 by an area bounded by two S curves. One has then 
to compare two broadly defined S areas rather than 
two sharply defined S curves. This can make an 
overlap of ordinates possible near the maximum of 
the (12, 6) curve and hence make solution of (3) possible 
where it previously was not. Furthermore, in com- 
paring the two S areas one has the possibility of choos- 
ing the two S curves, one within each area, whose 
slopes are most nearly alike. This could result in 
a reduction in the rapid variation of 62/61 with T* 
at low temperatures. For application to inherently 
imprecise experimental data, the low temperature 
region therefore becomes a much less sensitive probe 
of the potential than is indicated in figure 5. That 
is, an approximate fit to an accuracy compatible with 
experimental error might be possible where an exact 
fit, as indicated by figures 1 and 5, is impossible or, 
at best, difficult. 

There is another, more fundamental reason why the 
rapid variation of €2/61 with T* at low temperatures 
does not necessarily mean a sensitivity to differences 
in the potential functions. In this calculation, we 
have required the equality of both the second virial 
coefficient and its slope for the two potentials. This 
applies a much more stringent condition on the func- 
tions than is required in the correlation of experi- 



mental data. In the latter case, it is asked only that 
the theoretical values of B(T) come as close as possible 
to the experimental ones. Nothing is asked of the 
slope of the second virial coefficient. Clearly two 
functions may each fit the data within experimental 
precision, yet their slopes may disagree by con- 
siderably more than the precision of the present 
calculation. 

On the other hand, the fact that we place such strong 
conditions on the potential enables us to make strong 
statements where the second virials and their first 
derivatives for the potentials are essentially indis- 
tinguishable from each other. Obviously where our 
calculations cannot distinguish between potentials, 
a correlation which makes use of experimental data 
will be able to distinguish between them to a much 
lesser extent. Clearly, therefore, experiments de- 
signed to measure the second virial coefficient for 
purposes of learning something about the potential 
function should never be carried out above 7"* = 2.0 
on the (12,6) scale. In fact, existing data in that range 
should not be included in a determination of potential 
parameters since such data will supply experimental 
error without supplying any discrimination and so 
will reduce, for example, the ratios of the standard 
deviations obtained from fits of different potentials. 
This reduced temperature is easily converted to real 
temperatures for a particular substance given the 
ejk value for the (12,6) potential for that substance. 
For example, for argon, the data must have been taken 
at r<240°K while for xenon, the corresponding re- 
quirement is T < 450 °K. 

It is clear from figures 5 and 6 that in each class 
(i.e., square well, exp-6, etc.) there exists a potential 
for which the ratio 62/61 is essentially independent of 
7 1 * even at low temperatures. For the exp-6 this oc- 
curs for a slightly larger than 13. For the square well, 
it occurs for /?* approximately equal to 1.82. One 
expects this also to be true for other classes of three 
parameter potential classes of which the (12,6) is not 
a member. That is, there will exist a member of each 
such class which is equivalent to the (12,6) in pre- 
dicting the second virial coefficient over a large tem- 
perature range including low temperatures. Since 
the (12,6) potential function was chosen as the refer- 
ence potential in an entirely arbitrary fashion, there 
is no need to restrict this result to it. Thus, one can 
actually state that given any potential function, it is 
possible to find in every three parameter family of 
functions of which it is not a member, a potential 
function with which one can obtain a classical second 
virial coefficient whose value and slope differ from 
those calculated with the given potential by an amount 
much less than the best available experimental pre- 
cision over a temperature range starting at extremely 
low temperatures and extending to temperatures 
well above experimental conditions for almost all 
substances. In short, the second virial coefficient 
is seen to be at best a three parameter quantity with 
regard to the potential function and any attempt to 
use functions with more parameters necessarily leads 
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to redundancies. This is presumably what is behind 
the inability to obtain unique parameters in recent 
attempts to determine the potential function from 
second virial coefficient data using many parameter 
potential functions. 

As expected, deviations from this can occur at high 
temperatures. That is, where the repulsive parts 
of the potentials differ sufficiently in character, the 
high temperature region can be used to choose among 
different classes. Thus a choice can be made between 
the extremely different (12,6) and square well poten- 
tials if the data covers a range above T* = 7.0 on the 
(12,6) temperature scale. A choice between the some- 
what less different exp-6 and (12,6) potentials, on the 
other hand, requires data at T* > 10.0 on the same 
scale. These are very high temperatures for most 
substances. For argon this latter requirement is 
T> 1200 °K while for xenon it is T> 2250 °K. 

The ratio (60)2/(^0)1 for the second virial coeffi- 
cient behaves in essentially the same way as does 
e 2 /ei. There is, therefore, no need to discuss its 
behavior separately. 

Attempts to select, from several functions, a po- 
tential function for a particular system have sometimes 
been based on fits to second virial coefficient data 
which lie almost entirely within the flat portions of 
figures 5 to 8. Examples are the correlations of 
Whalley and Schneider [13] and of Mason and Rice 
[14]. In table 1 we have reproduced the standard devi- 
ations obtained by Whalley and Schneider for several 
potential functions. Note that only in the case of 
xenon, where half of the points lie outside the flat 
portions of figures 5 to 8, is there a strong discrimina- 
tion among the potential functions. For the krypton 
data, one certainly has no basis for the selection of 
one potential over another while for argon the choice 
is, at best, a marginal one. 



Table 


1. Standard deviations obtained by Whalley and Schneider a 
















Total 


Number of 




(9.6) 


(12,6) 


a=12 


a=13 


a=14 


a =15 


number of 
isotherms 


isotherms 

r < 2.0 

(12,6) 


Argon 


0.39 


0.42 








0.53 


15 


2 


Krypton 


.35 


.36 


0.46 


0.47 


0.41 


.42 


9 


2 


Xenon 


1.20 


.82 


.88 




.70 


.62 


12 


6 



a E. Whalley and W. G. Schneider, J. Chem. Phys. 23, 1644 (1955). 

From our results one can also see the futility of 
basing the choice of a potential function on the basis 
of the best fit of experimental data to a single two 
parameter function. To demonstrate this, let us 
take as an experimental system that system whose 
intermolecular potential function is exactly the (12,6) 
function. Figure 6 then represents an attempt to 
fit the "experimental" second virial coefficients to 
those predicted for the exp-6 potential. The best 
fit is obtained for that potential which gives the most 
nearly flat curve in figure 6. According to that figure, 
this best fit occurs for a value of a slightly greater 
than 13. The potential function defined by that value 



of a and the pair of parameters which give this best 
fit can then be associated with our "experimental" 
system and possibly used as such in other theories. 
Suppose now that instead of doing fits for a series of 
values of a we had just done the fit for a single value 
of a. Clearly, for every value of a, a pair of param- 
eters exists which gives the best fit to the "experi- 
mental" data for that value of a. However, the 
potential represented by that value of a and this 
pair of parameters could not in general, be associated 
with the experimental system unless the application 
is to a theory only weakly dependent on the potential, 
since the second virial coefficient associated with 
that function does not properly represent the low 
temperature second virial coefficient data. It is clear 
from figure 6, therefore, that one must take the best 
fit of experimental data to a series of two parameter 
potentials (here the family of functions generated 
by varying the third parameter, a) before assigning 
a particular potential to the experimental system. 
Unfortunately, the literature is full of fits of data to 
single potential functions, particularly to the (12,6) 
potential. Quite often the resulting potential has been 
used as the intermolecular potential functions for 
particular systems in evaluations of theories [15]. 
According to our results one must be suspicious 
of conclusions as to the relation between the par- 
ticular theory and experiment, based on such work, 
unless the theory is known to depend only weakly 
on the potential function. 



4. Zero Density Adiabatic Joule-Thomson 
Coefficient 

This quantity is obtained directly as the zero density 
limit of experimental free expansion data. To ob- 
tain the second virial coefficient, on the other hand, 
one must first, in some manner, differentiate the ex- 
perimental P-V-T data with respect to the density 
and then take the zero density limit. As a result, 
given the same experimental precision, one obtains 
the latter with much less precision than the former. 
This is not a real advantage for the free expansion 
data at the present time, however, since such data 
can be obtained only with a precision orders of mag- 
nitude below that possible in P-Y-T work, particu- 
larly at low densities. 

The adiabatic Joule-Thomson coefficient depends 
both on the second virial coefficient and its first de- 
rivative. Clearly, where two second virial coeffi- 
cients, one for each of two potentials, are indistinguish- 
able over an extended temperature range, their first 
derivatives are also indistinguishable, at least for 
temperatures near the center of this range. Further- 
more, this will be true for exactly the same parameter 
ratios. Differences which occur for the second virial 
coefficients near the edge of this temperature range 
necessarily appear as larger differences in their de- 
rivatives. Thus, one expects the range of equivalence 
for two potentials to be smaller for the zero density 
Joule-Thomson coefficient than it is for the second 
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virial coefficient. Results for this quantity are con- 
tained in figures 7 and 8. As expected, the range of 
equivalence of the potentials is shorter. It should 
be remembered that, since our method involves 
equating a property and its first derivative, equiva- 
lence here includes the second derivative of the second 
virial coefficient. 

5. Zero Density Viscosity and Diffgsion 
Coefficients 

The 62/61 ratios for the zero density viscosity and 
diffusion coefficients are contained in figures 6, 7, 9, 
and 10. The ratios (60)2/(60)1 for several potentials 
are presented in figure 11. In figure 7, for the Kihara 
potential, there is no curve for either of these proper- 
ties which approaches that of the second virial coef- 
ficient in flatness. The ratio €2/61, can be called 
independent of T* only in a very narrow temperature 
region about the maximum of the curve. Figure 9 
(which is drawn to a different scale) does indicate 
an increasing degree of flatness with increasing y. 

KIHARA POTENTIAL FUNCTION 



- 


1 1 1 

X--0.0 


1 1 1 1 


- 




~~~~~~~~~-~^_^ = ^ ::: ~ ::r ^ ::: ~ :: ^^---~^ 


/^ Z ^ EE1 




x=o.i 1 / ^~~ 


^^^Hr^^ 


0.2 / j / ' y^ 

°A/ L /-- 

N ^0.5 y 

VISCOSITY 

DIFFUSION 

1 1 1 


^^— -^ 


^^^^^-^ 


0.6 

0.6 

1 1 1 1 



3.0 4.0 5.0 

T*(I2,6) 



FIGURE 9. The parameter ratios 62/61 for the diffusion and vis- 
cosity coefficients for the Kihara potential function with respect 
to the (12,6). 

Notice the increasing flatness with increasing y for T* large. 
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Figure 10. The parameter ratios € 2 /ei for the diffusion and vis- 
cosity coefficients for the exp-6 potential with respect to the (12,6) 
potential. 



FIGURE 11. The parameter ratios (boW(bo)i for the Kihara and 
exp-6 potentials with respect to the (12,6). 



This occurs for relatively high temperature, T> Tr, 
however. This demonstrates the strong role played 
by the bowl of the potential function in determining 
the transport properties. For any given potential 
function, there is a temperature above which the 
transport properties are entirely dependent on the 
repulsive part of the potential function. Now, for 
the Kihara potential function, the position of the hard 
sphere cutoff moves to larger values of r as y in- 
creases. Therefore, the temperature at which the 
properties for the Kihara potential become those for 
a hard sphere should decrease with increasing y. 
According to figure 9, this is compensated for by the 
changes produced by the modification of the bowl with 
increasing y. The net result is a set of properties 
(e.g., for y = 0.6) more nearly like those for the y = 0, 
(i.e., (12,6)) potential at the highest temperatures of 
figure 9. Further proof of the strong effect produced 
by the bowl can be seen in the behavior at low tempera- 
tures. At such temperatures, the repulsive part has 
a very minor effect on the properties. Therefore, 
the rapid variation in the dependence of 62/61, with 
7 1 * as a function of y at low temperatures is an indica- 
tion of the marked effect of the bowl on these 
properties. 

The curves in figure 10 are most flat for T* approxi- 
mately equal to 2.0. It is interesting to note that 
there is a value of a for which the curve is quite flat 
at low temperatures as well. This occurs for a 
slightly larger than 13 for both properties. It should 
be noted that the corresponding ratio is very close 
to unity. Furthermore, for essentially the same value 
of a, a flat curve with essentially the same ordinate 
obtained down to low temperatures for the second 
virial coefficient. Note that, unlike the case of the 
second virial coefficient, for these properties the 
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curve for that value of a is not flat to very high tem- 
peratures. For these properties, a significant depar- 
ture from unity occurs at !T* = 5.0, not an impossibly 
high temperature for many systems. There is a 
tendency for the curves to flatten out at high tempera- 
tures for a value of a larger than 15 under which con- 
ditions the curves show a marked deviation from 
flatness at quite low temperatures. This behavior 
shows that it may be possible to distinguish between 
the (12,6) potential and all members of the family 
of exp-6 functions given experimental transport 
data of sufficient accuracy which includes both the 
temperature ranges T 7 * < 2.0 and 7 1 * > 5.0 on the 
(12,6) scale. For argon this requires that there be 
data for T< 240 °K and T > 600 °K while for xenon 
T<450 °K and T> 1225 ° K. These conditions are 
met for argon but not for xenon, at the present time. 

It is clear from these results that the diffusion and 
viscosity coefficients are potentially more sensitive 
probes of the potential function than either the second 
virial or Joule-Thomson coefficients given experi- 
mental data covering a sufficiently large temperature 
range. One might not have expected this since the 
latter properties are more simply related to the po- 
tential function than are the former. One might 
have expected a quantity like the second virial co- 
efficient which is related to the potential function 
through a single integration to be much more sensi- 
tive to differences in that function than one like the 
transport properties which are connected by three 
integrations, particularly when the functional depend- 
ence in the integrand is also more indirect for the 
latter than it is for the former. 

It must be remembered that the requirement here 
has been that both the value and slope of the prop- 
erties be equal for the two potentials. It is possible 
that some of the apparent sensitivity found for the 
transport properties comes from the requirement 
on the slope. Whether or not this disappears when 
one asks only for a match to the property, as is done 
in a correlation, bears further study. Again one can 
state that a lack of uniqueness found in this calcula- 
tion will not disappear on the application to experi- 
mental data. Therefore, one must have data present 
in the temperature ranges mentioned if one is even 
to have a chance to discriminate between potential 
functions. 

The slope of the €2/61 versus 7 1 * curve, for a given 
potential, is essentially the same for both the viscosity 
and diffusion coefficients. On this basis, the two quan- 
tities are equally effective when used individually. 
As we shall see in the next section, however, there 
is a difference between them in sensitivity to changes 
in the potential function. That difference makes 
the simultaneous fit of the two properties potentially 
a sensitive tool for finding the potential function. 

The ratios (60)2/(60)1 for these properties are less 
dependent on temperature than are the ratios e-z/ei. 
A discussion of them would therefore not contribute 
anything new with regard to the sensitivity of these 
properties as probes of the potential function. 



6. The Simultaneous Fit to More Than One 
Property 

Figures 6, 7, 9, and 10 contain plots of the ratio 
€2/61 for more than one property. With the help of 
these one can examine the effect the requirement 
of a simultaneous fit to two or more properties and 
their first derivatives has on the lack of uniqueness 
in the potential function. Figure 9 is particularly 
informative in this regard. Note that in that dia- 
gram, the curves for both the viscosity and diffusion 
coefficients for y=0 would be straight lines parallel 
to the abscissa at 62/61 = 1.0. As 7 deviates from zero, 
however, each curve begins to deviate from a straight 
line, particularly at low temperatures. Furthermore, 
the curves for the diffusion coefficient separate from 
those for the viscosity until, by 7 = 0.6, the curves 
are quite widely separated. The lack of uniqueness 
associated with the flatness of the curves for the 
separate properties would be extended to the simul- 
taneous fit to the pair of properties and their deriva- 
tives were the two curves essentially flat at the same 
value of the ordinate. However, the increase in the 
separation of the curves with y (i.e., with increased 
deviation from the (12,6) reference potential) dem- 
onstrates a strong sensitivity of the simultaneous 
fit to changes in the potential function. This sen- 
sitivity is much stronger than that of the individual 
properties in the range 7 1 * > 4.0. Whether or not 
it is more sensitive than the behavior for T* < 4.0 
would require an examination of detailed fits to data. 

Figure 10 contains the same kind of information 
for the exp-6 function. Here the changes that occur 
at low and intermediate temperatures are much smaller 
while those at high temperatures are much more 
pronounced. There appears to be a value of a for 
which the curves both have a flat range and are es- 
sentially superimposed at high temperatures. How- 
ever, for that value of a, the low temperature data 
still serves as a discriminant. In this case the 
simultaneous use of both properties and their deriva- 
tives does not appear to add any discrimination be- 
yond that available with either property by itself. 
Note that the curves for the two properties exchange 
relative positions when one goes from a = 12 to a = 17. 
Increasing a beyond 17 would presumably result in 
further separation. 

Figure 7 contains plots for all of the properties for 
the Kihara potential. These curves necessarily each 
form a straight line parallel to the abscissa with ordi- 
nate 1.0 for y = 0. Increasing y causes them to 
separate as in figure 9 for the two transport properties. 
It is interesting to note that the change in separation 
between the curves for diffusion and second virial 
coefficients at intermediate temperatures is not uni- 
form. Thus, for T* = 3.0, the diffusion curve for 
y = 0.1 lies above that for the second virial coefficient, 
for 7 = 0.3 they are superimposed in a small region 
around 7"* = 3.0, while for 7 = 0.4 the diffusion curve 
lies below. Thus, as 7 is increased from 0. to 0.1, 
the second virial curve "moves" more rapidly while 
for 7 > 0.1 the reverse is true. 
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7. Conclusions 

We have investigated the use of the second virial 
coefficient as well as the zero density Joule-Thomson, 
viscosity, and diffusion coefficients as probes for the 
intermolecular potential function. We have found 
| the second virial coefficient to be particularly poor 
in this regard in the temperature range 2.0 < 7 1 * < 7.0 
on the Lennard-Jones (12,6) reduced temperature 
scale. The Joule-Thomson coefficient has been 
found to be somewhat better particularly in the lower 
part of this range. Considerably better were the 
transport coefficients, i.e., the viscosity and diffusion 
coefficients. It is clear from our results that these 
last coefficients can be sensitive probes of the poten- 
tial function given experimental data covering a suf- 
ficiently large temperature range. 

Further comparisons of this kind between the 
transport coefficients for the (12,6) potential function 
and those for other potential functions are clearly 
called for, particularly for classes of potentials of 
which the (12,6) is not a member. In particular, 
one would like to find such classes of functions as 
exist which contain a member whose transport prop- 
erties can be fit to those of the (12,6) over an extensive 
temperature range. This will require the calculation 
of collision integrals for these various functions. 
In some cases tables do exist [16, 17J which, unfortu- 
nately, contain too few points in the temperature 
region of interest. It is hoped that these will be ex- 
tended shortly. 

These calculations will be extended to include the 
thermal diffusion ratio and, possibly, the quantum 
corrections to certain of these properties. 

As regards the relationship between our results 
and experiment with present day precision, the fact 
that our results apply precisely to experiment only 
in the limit of very high precision actually strengthens 
our conclusions in regard to lack of uniqueness while 
conclusions regarding uniqueness are weakened. 
That is, where we did not find it possible to use a 
macroscopic property to distinguish between poten- 
tial functions under our conditions, it certainly would 
not be possible to use this property for this purpose 
under less precise experimental conditions. On the 
other hand, our ability in other circumstances to dis- 
tinguish among potentials (e.g., using second virial 
coefficient data for T* (12,6, < 2.0) may be due in part, 
to the fact that we have required a precise fit, some- 
thing not possible with experimental data. It may 
likewise be due to the inclusion here of a requirement 
that the first derivative of the property be equal for 



the two potentials. More precise statements in this 
latter case await the result of fits to actual data. 
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